International  Journal  of  Power  Electronics  and  Drive  System  (IJPEDS) 

Vol.  8,  No.  4,  December  2017,  pp.  1663-1672 

ISSN:  2088-8694,  DOI:  10.11591/ijpeds.v8i4.ppl663-1672 


□  1663 


Permanent  Magnet  Tubular  Generator  with  Quasi-Halbach 
Array  for  Free-Piston  Generator  System 


Behrooz  Rezaeealam 

Departement  of  Electrical  Engineering,  Lorestan  University 


Article  Info 


ABSTRACT 


Article  history: 

Received  Sep  1,  2017 
Revised  Sep  22,  2017 
Accepted  Nov  3,  2017 


Keyword: 

PM  linear  generator 
Finite  element  method 
Imperialist  competitive 
algorithm 


In  this  paper,  the  performance  of  a  free-piston  generator  system  with  tubular 
single-phase  stator  winding  is  investigated  for  three  different  configurations 
of  permanent  magnet  arrays  on  the  mover:  radial,  rectangular  quasi-Halbach 
and  trapezoidal  quasi-Halbach  arrays.  The  optimal  shape  of  the  rectangular 
quasi-Halbach  magnetization  is  obtained  by  employing  an  optimization 
technique  called  Imperialist  Competitive  Algorithm  (ICA),  and  the  results 
are  compared  with  the  previously  developed  permanent  magnet  linear 
generator  (PMLG)  with  magnets  of  radial  magnetization.  Then,  the  resultant 
design  of  rectangular  quasi-Halbach  array,  is  changed  to  trapezoidal  quasi- 
Halbach  array  and  its  performance  is  compared  with  other  types  of  magnet 
arrays.  The  employed  ICA  is  based  on  the  analytical  model  of  the  magnetic 
flux  distribution  in  air-gap,  and  is  further  verified  by  a  nonlinear  transient 
finite  element  (FE)  model.  The  FE  model  has  been  developed  to  study  the 
performance  of  the  reciprocating  PMLG.  The  accuracy  of  the  proposed  FE 
model  is  exemplified  by  comparing  the  calculated  results  with  the 
experimental  ones  regarding  the  PMLG  with  radial  magnetization. 
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1.  INTRODUCTION 

High  oil  prices  have  put  the  environmental  friendly  cars  in  the  spotlight.  Now,  many  car 
manufacturers  are  developing  some  alternative  solutions  to  the  conventional  car.  A  free-piston  generator  is 
employed  for  energy  conversion  in  which  an  internal  combustion  engine  and  a  linear  electric  generator  are 
integrated  into  a  single  unit  [1,2].  Hence,  the  crankshaft  and  the  connecting  rod  are  omitted.  As  the  moving 
parts  of  the  engine  are  minimized,  the  flexibility  and  efficiency  of  the  generator  are  improved.  The  free- 
piston  generator  is  a  proper  choice  for  the  series  hybrid  electric  vehicles,  because  the  fuel  consumption  is 
lower,  and  thereby,  the  environmental  pollution  may  decrease.  Also,  the  interesting  feature  of  the  free-piston 
generator  is  the  multi-fuel  functioning  capability  [3]. 

As  the  efficiency  and  power  density  of  PMLG  are  higher  than  the  ones  of  induction -type  and 
switched  reluctance  generators  [4, 5, 6, 7],  in  recent  years,  several  types  of  PMLG  have  been  developed  for 
free-piston  generator  system.  The  PMLG  with  transverse  flux  topology  has  high  magnetic  flux  density, 
however,  it  suffers  from  high  leakage  flux  [6] .  The  PMLG  with  radial  magnetization  that  has  been  employed 
for  free  piston  generator  system  [1,2],  has  a  detent  force  due  to  the  interaction  of  the  mover  permanent 
magnets  with  the  stator  core.  Generally,  the  intense  detent  force  adds  to  the  thrust  ripple  of  the  PMLG. 

The  PMLG  with  quasi-Halbach  topology  has  a  self-shielding  property  and  higher  flux  density  than 
the  one  with  radial  topology  [8,9].  Furthermore,  by  replacing  the  rectangular  shaped  magnets  of  quasi- 
Halbach  array  with  trapezoidal  shaped  magnets,  the  average  flux  density  in  air  gap  of  tubular  actuators  is 
elevated  [10].  In  literature,  optimal  design  of  quasi-Halbach  topology  has  so  far  been  presented  based  on 
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different  modeling  techniques  including  magnetic  equivalent  circuit  with  lumped  elements,  analytical  method 
using  Maxwell  equations,  and  finite -element  method  [11,12]. 

In  this  paper,  in  order  to  increase  the  thrust  force  and  output  power  of  the  free  piston  generator 
system,  the  optimized  PMLG  with  quasi-Halbach  magnetization  is  employed.  At  first,  optimization  of  the 
PMLG  is  performed  using  ICA,  that  is  based  on  the  analytical  model  of  the  PMLG,  and  an  effective  objective 
function  is  suggested,  regarding  the  thrust  force  and  the  flux  density  distribution  in  air  gap.  Next,  the 
performance  of  the  optimized  PMLG  with  quasi-Halbach  magnetization,  is  investigated  using  FE  method. 
For  this  purpose,  a  2-D  nonlinear  transient  FE  model  is  developed  to  consider  the  reciprocating  motion  of  the 
mover,  by  moving  mesh  techniques.  Finally,  it  is  shown  that  further  improvement  in  the  performance  of  the 
generator  can  be  achieved  by  replacing  the  rectangular  magnets  of  the  quasi-Halbach  topology  with  the 
trapezoidal  magnets. 


Figure  1.  Free  piston  generator  system  with  a  PMLG  [1]. 


2.  FREE-PISTON  GENERATOR  SYSTEM 

The  employed  combustion  engine  in  the  system,  is  a  linear,  crankless,  two  stroke  engine,  as  shown 
in  Figure  1  in  which  two  horizontally  opposed  pistons  are  fixed  on  a  common  connecting  rod  that  oscillates 
between  the  two  end-mounted  cylinders.  Ignition  happens  alternately  in  each  cylinder,  driving  the  piston 
back  and  forth.  The  bore  of  the  engine  is  36.5  mm  and  the  maximum  stroke  length  is  50  mm. 

The  profile  of  the  reciprocating  motion  of  the  mover  is  numerically  modeled  using  the  following 
force  balance  equation  [1]: 
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where  AB  ,  Pi  ,  r,  x,  xm,  n  and  m  are  bore  area,  intake  pressure,  compression  ratio,  mover  position  maximum 
half  stroke,  ratio  of  specific  heat  and  mass  of  the  mover,  respectively.  The  mass  of  the  mover  includes  the 
mass  of  the  moving  parts  of  the  PMLG. 

The  PMLG  consists  of  two  main  parts,  a  stator  and  a  mover.  The  stator  consists  of  the  armature  core 
and  the  single-phase  windings,  and  the  mover  is  made  up  of  the  PM’s  with  radial  magnetization  and  the  shaft 
on  which  they  are  mounted  as  shown  in  Figure  1.  When  the  engine  is  running,  the  mover  of  the  PMLG 
reciprocates  relative  to  the  stator  windings,  thereby  the  distribution  of  the  magnetic  field  in  air  gap  changes, 
leading  to  voltage  induction  in  the  stator  winding.  Thus,  the  output  power  could  be  extracted  from  the  PMLG 
as  the  load  is  connected  to  the  stator  windings.  The  geometric  parameters  of  the  PMLG  with  radial 
magnetization,  are  listed  in  Table  1. 


3.  PMLG  WITH  QUASI-HALBACH  MAGNETIZATION 

Figure  2  shows  the  PMLG  with  quasi-Halbach  magnetization,  besides  the  one  with  radial 
magnetization.  In  this  paper,  the  radial  topology  is  replaced  with  the  quasi-Halbach  topology  for 
reciprocating  free-piston  generator  system. 

The  distribution  of  the  radial  component  of  the  flux  density  Br  in  the  air  gap  along  the  axial  distance 
z  for  PMLG  with  quasi-Halbach  magnetization  can  be  expressed  as  [9,13]: 

Br  =  Xn=i,2,SamBli  (™nr)  +  b^BK-^  ( mnr ))  sin  (mnz)  (2) 
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where,  mn  =  (2 n  —  l)n/rp  and  Tp  is  the  magnet  pole  pitch,  B f  is  the  modified  Bessel  function  of  the  first 
kind  of  order  1  and  BK±  is  the  modified  Bessel  function  of  the  second  kind  of  order  1.  a/n,  bIn ,  a//n,  and  bIIn 
have  been  described  in  [13]  that  also  depend  on  the  geometry  of  the  PMLG. 

Furthermore,  the  thrust  force  resulting  from  the  interaction  between  the  winding  current  and  the 
permanent  magnet  field,  is  given  by: 

Fw  =  £”=1,2,...  Fn  sin  mnz  (3) 


Table  1.  PMLG  parameters 


Number  of  active  rotor  poles 

P 

4 

Poles 

Pole  pitch 

T 

41.91 

mm 

Diameter 

D 

62.865 

mm 

Number  of  windings 

Ns 

296 

turns 

Magnet  Thickness 

K 

18.415 

mm 

Magnet  width 

h 

41.0718 

mm 

Magnet  Gap 

wmg 

0.4191 

mm 

Tooth  Width 

Wt 

20.955 

mm 

Slot  Width 

ws 

20.955 

mm 

Back  Iron  Thickness 

hi 

14.7066 

mm 

Air  Gap 

9 

1.651 

mm 

Magnet  Inner  Diameter 

IDm 

26.035 

mm 

Magnet  outer  Diameter 

0Dm 

62.865 

mm 

Mover 

(secondary) 


Stator 

(primary) 


Shaft 


► 


1 


1 


(a) 


Mover 

(secondary) 


Stator 

(primary) 


(b) 


Figure  2.  PMLG  topologies,  (a)  Radial  array,  (b)  quasi-Halbach  array, 
where,  the  harmonic  components  Fn  have  been  described  in  [13],  and  the  effective  force  is  defined  as: 


Fnrms  =  JZTJt  (4) 

4.  OPTIMIZATION  OF  THE  PMLG  WITH  QUASI-HALBACH  MAGNETIZATION  USING  ICA 

The  PMLG  with  radial  magnetization  in  [1,2]  is  considered  as  the  basis  for  the  design  of  the  PMLG 
with  quasi-Halbach  magnetization.  As  depicted  in  Figure  2,  the  design  variables  are  the  height  of  magnet 
(/im),  the  width  of  the  magnets  (rmx,  rmy),  the  stator  slot  pitch  (rsp),  width  and  height  of  the  stator  slot 
(, b0 ,  hs),  and  the  number  of  series  turns  in  one  slot  ( Nsp )  in  the  PMLG  with  quasi-Halbach  magnetization.  The 
fixed  variables  are  pole  pitch  rp,  pole  pairs,  air  gap  length,  and  also,  the  length  and  the  outer  diameter  of  the 
generator  due  to  the  mechanical  constraints  in  the  mover.  The  pole  pitch  is  considered  fixed,  because  it  is  set 
equal  to  the  stroke  length  of  the  engine  to  reduce  the  harmonic  content  of  the  output  waveforms,  and  then: 
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Now,  the  optimization  technique  of  ICA  is  employed  to  obtain  the  optimal  shape  of  the  PMLG  with 
quasi-Halbach  magnetization,  by  maximizing  the  flux  density  in  the  air  gap  and  also  the  effective  thrust 
force.  For  this  purpose,  the  fitness  function  is  defined  as  follows  [14]: 

FF  =  (a1.Fnrms  +  a2.Brl)(xi . ^  (6) 

where  xlf  ...,xn  are  design  variables,  and  +  a2  =  1,  that  the  index  weights  of  0.6  and  0.4  are  considered 
for  a1  and  a2 ,  respectively. 

The  ICA  is  employed  to  search  for  the  maximum  value  of  FF.  The  ICA  provides  a  random  search 
method  to  get  the  overall  optimal  solution  in  a  complex  multi-dimensional  search  space.  In  this  algorithm,  the 
number  of  population  is  equal  to  the  number  of  countries.  Every  country  has  an  array  of  variables  that  is 
defined  as  follows: 

country  =  [hm,  rmx,  Tsp,  b0,  hs>  Nsp]  (7) 

In  order  to  start  the  algorithm,  some  of  this  population  should  be  created.  Thus  the  total  population 
matrix  is  created  randomly.  The  total  population  matrix  is  defined  as  follows: 


country1 

country2 


COUNTRY = 


[country  Ncountry\ 

^■mx1 

h"m2  ^mx2 
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K 


h-s2 
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°oNi 
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The  costs  of  the  countries  are  defined  as  follows: 


costi  =  f  (country  i)  =  f(hm,  rmx,  rsp,  b0,  hs,  Nsp) 


(8) 


(9) 


The  objective  function  is  defined  as  (6)  for  increasing  the  thrust  force  and  the  flux  density.  As 
illustrated  in  Figure  3,  the  ICA  with  parameters  listed  in  Table  2,  is  employed  to  search  for  the  optimal 
design.  The  optimal  design  is  described  by  the  resultant  values  for  the  parameters  listed  in  Table  3. 


Table  2.  ICA  parameters. 

symbol 

Quantity 

value 

Up  0p 

Number  of  countries 

500 

W-imp 

Number  of  imperials 

10 

TLVar 

Number  of  unknown  variables 

4 

d-dec 

Number  of  decades 

500 

Table  3.  Specification  of  optimal  designed  generator. 


parameter 

optimal 

hm 

18.415  [mm] 

^mx 

16.428  [mm] 

T  my 

24.64  [mm] 

hs 

14.706  [mm] 

^0 

20.955  [mm] 

Tsp 

41.91  [mm] 

Nsp 

39  [turns] 
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5.  TRANSIENT  FINITE  ELEMENT  MODELING 

In  this  section,  the  PMLG  with  radial  array  described  in  [1,2]  and  also,  the  one  with  quasi-Halbach 
array  structure  are  simulated  using  time-stepping  transient  FEM.  The  Multiphysics  finite  element  package 
COMSOL  considers  the  relative  movement  between  the  translator  and  the  stator  using  moving  mesh 
techniques  [15].  The  stator  mesh  remains  fixed  and  is  connected  to  the  moving  translator  mesh.  The  shapes 
of  the  generator  changes  during  the  solving  process  and  the  new  mesh  is  obtained  according  to  the  new  shape 
of  the  generator.  Figure  4  shows  the  triangular  mesh  for  initial  position  of  the  translator  and  also  at  the 
moment  when  the  translator  reaches  its  end  position. 


Figure  3.  ICA  flowchart. 


Initial  Mesh 


compressed  Mesh 


(a)  (b) 

Figure  4.  The  triangular  mesh:  (a)  initial  position  of  the  translator  (b)  the  moment  at  which  the  translator 

reaches  its  end  position 


The  finite  element  model  of  the  optimized  PMLG  with  quasi-Halbach  array  is  shown  in  Figure  5 
that  illustrates  the  distribution  of  magnetic  flux  density  inside  the  generator.  The  quasi-Halbach  topology  has 
a  self-shielding  property  and  higher  flux  density  in  the  air-gap  region  in  comparison  to  the  radial  topology  as 
shown  in  Figure  6.  Also,  the  distribution  of  magnetic  flux  density  is  more  uniform  inside  the  generator  with 
quasi-Halbach  topology. 
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Figure  5.  3D  view  of  the  magnetic  flux  density  distribution  in  PMLG 
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Figure  6.  Flux  density  distribution  in  the  generators:  (a)  Radial  array,  (b)  quasi-Halbach  array. 


The  design  optimization  in  this  work  was  carried  out  based  on  the  analytical  model  of  the  machine 
presented  in  Section  4.  As  the  validity  of  the  design  optimization  greatly  depends  on  the  accuracy  of  the 
model,  therefore,  for  the  purpose  of  comparison,  the  magnetic  flux  density  distribution  as  a  function  of 
position  along  the  air-gap  is  depicted  in  Figure  7  for  the  PMLG  with  quasi-Halbach  array,  that  has  been 
calculated  using  analytical  equation  of  (2)  and  also  using  the  finite  element  model  of  the  generator,  which 
demonstrates  good  agreement  between  them. 


Figure  7.  Comparison  between  analytically  predicted  and  FEM  calculated  distributions  of  the  flux 
density  in  air  gap  of  the  optimized  PMLG  with  quasi-Halbach  array. 
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Figure  8  shows  the  analytical  results  for  the  magnetic  flux  density  distribution  as  a  function  of 
position  along  the  air-gap  for  both  generators.  Higher  magnetic  flux  density  is  seen  inside  the  air-gap  of  the 
generator  with  quasi-Halbach  array  than  with  radial  array.  Figure  9  shows  the  thrust  force  for  radial  array  and 
quasi-Halbach  array  topologies  and  it  is  seen  that  trust  force  increases  effectively  in  a  PMFG  with  quasi- 
Halbach  array.  Figure  10  shows  the  output  power  of  the  generators.  For  the  case  of  PMFG  with  radial  array, 
the  derived  output  power  from  the  simulation,  matches  roughly  with  the  experimental  one  [1,2],  justifying  the 
suitability  of  the  developed  transient  finite  element  model.  The  engine-generator  system  was  capable  of 
supplying  a  peak  power  of  315W  for  radial  array  [1,2]  and,  over  the  same  range  of  load,  the  PMFG  with 
quasi-Halbach  array  provides  the  peak  power  of  336  W  as  shown  in  Figure  10. 


Figure  8.  Magnetic  flux  density  distribution  in  air  Figure  9.  Thrust  force  waveforms  of  the  radial  array 
gap  of  the  PMFG  with  quasi-Halbach  and  radial  and  the  quasi-Halbach  array  topologies 

array  topologies 


Figure  10.  Output  power  of  the  PMFG  with  radial  and  quasi-Halbach  arrays  versus  load  conductance 


6.  TRAPEZOIDAL  MAGNETS 

In  the  previous  section,  an  analytical  model  was  employed  to  optimize  the  design  of  the  quasi- 
Halbach  topology  with  rectangular  magnets  and  also,  a  finite  element  model  was  developed  to  examine  the 
performance  of  the  generator  in  more  detail.  In  this  section,  the  finite  element  model  is  extended  to  be  able  to 
investigate  the  influence  of  the  angle  am  of  the  trapezoidal  magnet  on  the  performance  of  the  PMLG,  as 
shown  in  Figure  1 1 . 

The  flux  density  distribution  in  air  gap  is  calculated  for  several  values  of  the  parameter  am  as 
indicated  in  Figure  12.  It  is  found  that  the  maximum  amplitude  of  flux  density  inside  the  air-gap  is  obtained 
at  am  =  110°,  however,  the  highest  average  value  of  the  flux  density  is  achieved  at  am  =  100°. 
Furthermore,  Figure  13  shows  the  output  power  of  PMLG  with  trapezoidal  magnets  for  different  values  of 
am  and,  it  is  noticed  that  the  maximum  output  power  of  359  W  is  attained  with,  am  =  100°. 

As  the  mover  of  the  PMLG  reciprocates,  the  distribution  of  the  magnetic  field  varies  along  the  axial 
direction  and  therefore  the  linkage  flux  of  the  stator  winding  varies,  leading  to  the  induced  back-EMF  in  the 
stator  winding.  Based  on  the  developed  finite  element  model,  the  back-EMF  waveforms  of  the  PMLG  are 
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shown  in  Figure  14,  for  the  three  cases  of  radial,  rectangular  quasi-Halbach  and  trapezoidal  quasi-Halbach 
( am  =  100°).  The  results  show  the  performance  improvement  with  trapezoidal  magnets. 


Figure  11.  Quasi-Halbach  magnetization  patterns: 
(a)  Rectangular  shaped  PM’s.  (b)  Trapezoidal 
shaped  PM’s 


Figure  12.  The  effect  of  am on  the  waveform  of  the 
radial  component  of  the  flux  density  in  air  gap 


Time  (s) 


Figure  13.  The  effect  of  amon  the  waveform  of  the 
output  power  of  the  PMLG  with  trapezoidal 
magnets 


Figure  14.  Back  EMF  waveform  of  the  three 
generators 


Figure  15.  The  translator  position  profile  for  both  Figure  16.  The  output  voltage  versus  position 

designs  of  radial  and  quasi-Halbach  arrays 
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Table  4.  Performance  of  radial  array  and  quasi-halbach  array  in  free  piston  generator. 


radial  array 

quasi-Halbach  array 
(Rectangular  shape) 

quasi-Halbach  array 
(Trapezoidal  shape) 

Rms  force  (N) 

2208 

2465 

2632 

Voltage  (v) 

66 

74 

85 

Load  power  (w) 

315 

336 

359 

There  is  a  gap  between  the  magnets  with  radial  array  as  illustrated  in  Figure  2,  while  the  magnets 
are  bonded  together  with  quasi-Halbach  array.  Because  of  this,  the  translator  with  quasi-Halbach  array  is 
heavier  than  the  one  with  radial  array.  Hence,  the  reciprocating  frequency  of  the  translator  with  quasi  - 
Halbach  array  (21.27  Hz)  is  slightly  lower  than  the  one  with  radial  array  (23.53  Hz)  [1].  The  profiles  of 
translator  position  are  shown  in  Figure  15  for  both  cases  of  radial  and  quasi-Halbach  arrays.  It  is  worth 
mentioning  that  the  translators  with  quasi -Halbach  array  with  either  rectangular  or  trapezoidal  shapes  have 
the  same  weight.  Also,  the  output  voltage  versus  position  is  shown  in  Figure  16  for  both  cases  of  radial  and 
quasi-Halbach  arrays.  The  increase  in  the  output  voltage  of  the  PMFG  is  seen  in  the  sequence:  radial, 
rectangular  quasi-Halbach,  trapezoidal  quasi-Halbach. 


7.  CONCLUSIONS 

In  this  paper,  PMLG  with  rectangular  and  trapezoidal  quasi-Halbach  magnet  arrays  were  proposed 
for  free-piston  generator  system.  The  ICA  was  employed  to  optimize  the  shape  of  the  rectangular  quasi - 
Halbach  array,  based  on  the  corresponding  analytical  model  for  the  air  gap  magnetic  flux  density  distribution 
and  thrust  force.  Also,  a  transient  FE  model  was  developed  by  which  the  performance  of  the  reciprocating 
PMLG  was  examined  by  incorporating  moving  mesh  techniques.  The  results  from  the  proposed  numerical 
model  of  the  PMLG  with  radial  array  were  in  good  agreement  with  experiments. 

In  addition,  the  rectangular  quasi-Halbach  array  was  changed  to  trapezoidal  quasi-Halbach  array, 
and  the  performances  of  the  three  mentioned  PMLG  were  compared  as  shown  in  Table  4.  The  results  show  a 
performance  enhancement  of  the  PMLG  with  quasi-Halbach  trapezoidal  magnets  compared  to  the  PMLG 
with  radial  magnetized  topology,  by  increasing  the  output  power  up  to  13%. 
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